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We analyze the enhancement of the superconducting critical temperature of superconduct- 
ing/ferromagnetic bilayers due to the appearance of localized superconducting states in the vicinity 
of magnetic domain walls in the ferromagnet. We consider the case when the main mechanism of 
the superconductivity destruction via the proximity effect is the exchange field. We demonstrate 
that the influence of the domain walls on the superconducting properties of the bilayer may be quite 
strong if the domain wall thickness is of the order of superconducting coherence length. 
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I. INTRODUCTION 

The coexistence of singlet superconductivity and fer- 
romagnetism is very improbable in bulk compounds but 
may be easily achieved in artificially fabricated hybrid 
superconductor (S)-ferromagnet (F) structures. 

There are two basic mechanisms responsible for inter- 
action of superconducting order parameter with magnetic 
moments in the ferromagnet: the electromagnetic mech- 
anism (interaction of Cooper pairs with magnetic field 
induced by magnetic moments) and the exchange inter- 
action of magnetic moments with electrons in Cooper 
pairs. The second mechanism enters into play due to the 
proximity effect, when the Cooper pairs penetrate into 
the F layer and induce superconductivity there. In S/F 
bilayers it is possible to study the interplay between su- 
perconductivity and magnetism in a controlled manner, 
since we can change the relative strength of two compet- 
ing orderings by varying the layer thicknesses and mag- 
netic content of F layers. Naturally, to observe the in- 
fluence of the ferromagnctism on the superconductivity, 
the thickness of S layer must be small. This influence 
is most pronounced if the S layer thickness is smaller 
than the superconducting coherence length £ s . Recently, 
the observation of many interesting effects in S/F sys- 
tems became possible due to the great progress in the 
preparation of high-quality hybrid F/S systems - see the 
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reviewsifii*. 

In practice, the domains appear in fcrromagnets and, 
near the domain walls, a special situation occurs for prox- 
imity effect. For the purely orbital (electromagnetic) 
mechanism of the superconductivity destruction, the nu- 
cleation of the superconductivity in the presence of do- 
main structure has been theoretically studied in^ for the 
case of magnetic film with perpendicular anisotropy. The 
conditions of the superconductivity appearance occur to 
be more favorable near the domain walls due to the par- 
tial compensation of the magnetic induction. Recently, 
the manifestation of such domain wall superconductiv- 
ity (DWS) was revealed on experiment where Nb film 
was deposited on top of the single crystal ferromagnetic 
BaFei20ig covered with a thin Si buffer layer. 

As the typical value of the exchange field in the fer- 



romagnets h ~ (100 — 1000)1^ exceeds many times the 
superconducting critical temperature T c o, the exchange 
mechanism prevails the orbital one in the superconduc- 
tivity destruction when the electrical contact between S 
and F layers is good. For the proximity effect mediated 
by the exchange interaction, the Cooper pairs feel the 
exchange field averaged over the superconducting coher- 
ence length. Naturally, it will be smaller near the domain 
wall and we may expect that superconductivity would be 
more robust near them. The local increase of the critical 
temperature in presence of magnetic domains was ob- 
served experimentally in Nio.soF e o.2o/Nb bilayers (with 
Nb thickness around 20 nm)l, and it was attributed to 
DWS formation. 

In the present paper we study theoretically the condi- 
tions of the localized superconductivity appearance near 
the domain wall taking into account exchange mecha- 
nism of the proximity effect. In Sec. [HI we demon- 
strate that, in the case of thin F layer and small domain 
wall thickness, the problem is somewhat similar to that 
of the domain wall superconductivity in ferromagnetic 
superconductor^. For the case when the superconduct- 
ing coherence length £ s exceeds the DW thickness w, we 
expect a very strong local increase of T c (see Sec. Illlll . 
In Sec. IIVI we obtain the analytical expression for the 
critical temperature of the DWS for the case when the 
DW thickness w exceeds the superconducting coherence 
length £ s . The Appendix presents an extension of this 
result to arbitrary thickness of the F layers and trans- 
parency of the S/F interface. We discuss our results in 
Sec. In particular, we predict the realization of the 
situation when superconductivity appears only near the 
DW. 



II. EQUATION FOR THE CRITICAL 
TEMPERATURE IN THIN BILAYERS 

We introduce the Usadel equations^ which are very 
convenient when dealing with S/F systems with critical 
temperatures T c and exchange fields h such as T c t <C 1 
and hr <C 1, where r is the elastic scattering time. 

Near the second order transition into superconducting 



2 



state, Usadel equations can be linearized with respect 
to the amplitude of the superconducting gap. In the S 
region, the linearized Usadel equation is: 



D s V 2 f s + 2\uj\f s = 2A s (r)a z , 



(1) 



where D s is the diffusion constant in the superconductor 
and ui = (2n + l)7i\T is a Matsubara frequency at tem- 
perature T. (Notations are similar to the ones used in 
Ref. 0, except for the factor i in front of A s ). Eq. 
relates the anomalous Green's function f s , which is a 
matrix in spin space, to the superconducting gap A s (a;), 
&(x.y,z) are Pauli matrices in spin space. In the absence 
of supercurrent, the gap can be taken as real. 

In F region, an exchange field hf — {hf^ x ,hf^ y ,hf^ z ) 
is acting on the spins of conduction electrons and the 
linearized Usadel equation for anomalous function ft is 

-D f V 2 f f + 2\w\f f + is w (hf, x {a x ,f f ] 

+h/,v[&v>ff] + hfA**Jf}) = °- (2) 

Here, Df is the diffusion constant in the ferromagnet and 
we used the abbreviation s w = sgn(w). 

In addition, the boundary conditions at the interfaces 
with vacuum yield d z f s (z = d s ) = and d z ff(z = 
—df) = 0, where d s and df are the thicknesses of S and F 
layers, respectively, and z = defines the plane of the in- 
terface between both layers. We also consider that S and 
F layers are separated by a thin insulating tunnel bar- 
rier. Therefore, the boundary conditions at the interface 
z = areiii 

<T s dJ s (0) = <TfdJf(0), f s (0) = //(0)+7 B & djf 

"(3) 

where a s and a f are conductivities in the layers, and 7s 
is related to the boundary resistance per unit area Rb 
through jbCs = Rb&f, where £ s is the superconducting 
coherence length. 

The critical temperature T — T c at the second order 
transition is now obtained from the self consistency equa- 
tion for the gap: 



Afrlln- 



T 



(•0 



A(r) 



-f^(r,u) =0 (4) 



(J 11 is a matrix element of /), where T c o is the bare 
transition temperature of the S layer. 

In the ferromagnet, the magnitude hf of the exchange 
field hf is fixed. However, its orientation can rotate in 
the presence of a magnetic domain wall structure. In 
the following, we assume a one-dimensional domain wall 
structure, along a;-axis. In order to find the critical tem- 
perature of the bilayer in the presence a domain wall, 
we must find the x- and z-dependence of the gap and 
anomalous functions f s and ff which solve Eqs. QJ- 
Proximity effect can significantly affect the tran- 
sition temperature only when the thickness of S layer is 



compar able with the superconducting coherence length 
£ s = y/ D s /2irT c o. In order to get tractable expressions, 
we will only consider the case d s <C £ s . Such regime 
is also well achievable experimentally. Then, f s and A s 
are almost constant along z-axis. Therefore, we can av- 
erage Eq. JIJ on the thickness of S layer and make use 
of the boundary condition at the interface with vacuum. 
Finally, we get the following equation at z = 0: 



D s d 2 j s 



d z f s +2\u\f s = 2A s a z 



(5) 



The characteristic scale for the proximity effect in F 
layer is rather set by coherence length £/ = y/Df /hf, 
where hf is the typical amplitude of the exchange field 
in the ferromagnet. In this section, we will address the 
case of very thin F layer: df <C £/. Then, from Eq. 
we can derive similarly the Usadel equation averaged over 
the thickness df, at z = 0: 

- Dfdlff - ^djf + 2\w\ff + is u (h i ,a,.f i (6) 



+hf, y [a y ,ff]+hf !Z {a z ,ff}J = 0. 

Let us note right now that the assumption of very thin 
F layer is quite hard to achieve experimentally, as we 
will discuss at the end of Sec. El I n Appendix, we will 
consider the case of arbitrary thickness for the F layer. 

For simplicity, we also consider the case of low interface 
resistance (7^ — * 0) where the proximity effect is max- 
imal. In this regime, ff{x) « f s {%) = f(x). By proper 
linear combination Eqs. 10 El EJ), we can form a single 
equation on f(x): 



- Dd 2 J + 2M/ + is u (h x [& x ,ff] + h y [a v , ff] 



+h z {a z Jf}) =2Aa z , A 



Vs + Vf 



A S! (7) 



where rj s = a s d s /D s and rjf = ufdfjDf. Therefore, 
the thin bilayer is described by the same equations as 
for a magnetic superconductor^ with effective diffusion 
constant, exchange field, and BCS coupling constant: 



D 



D s r) s + D f 7]f 
Vs+Vf 



h = ^^h f , A = — — — A, 

Vs + Vf Vs + Vf 

.(8) 

respectively. An equation similar to Eq. Q was derived 
for a thin normal- metal/superconductor bilay er, in the 
absence of exchange field (h = 0) in Ref. [l^. There, it 
was shown that the reduction of the coupling constant A 
leads to a rapid decrease of the bilayer critical tempera- 
ture. In the following, we do not consider these effects. 
Rather, we dwell with the case when A f=a A and the re- 
duction of T c is mainly due to effective exchange field h. 
Such situation occurs at r) s 3> ??/. (When S and F layers 
have comparable diffusion constant and conductivity, the 
renormalization factors in Eq. JHJl receive a simple inter- 
pretation in terms of volume ratios. In particular, the 



Vs 
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condition rj s ^S> T]f results in d s ^> df.) Thus. D w D s 
and h « {rjf /r) s )hf. Let us note that the amplitude of h 
is strongly reduced compared to ft./ , eventually it is of the 
order of A s , and thus it leads to the possible coexistence 
of magnetism and superconductivity in the bilayer. 

The phase diagram of magnetic superconductors with 
constant exchange held was studied long &goi& Second 
order transition line from normal to superconducting 
state at the critical temperature T = T c (h) is given by 
equation 



ln-^+27rTReV (- 



ih 



0, 



(9) 



where h is the amplitude of h. At zero temperature, the 

( 2) 

critical field is h c = A /2, where A ~ 1.76T c0 is the 
superconducting gap. However, at T < T* « 0.56T c0 , 
the transition into superconducting state is of the first 
order and the critical field at zero temperature is rather 
h c = A /V2. 

In the presence of a domain structure in the ferromag- 
net, the average exchange field felt by the electrons near 
domain walls is smaller than in the domains. This may 
lead to the enhancement of the superconducting critical 
temperature. On the basis of the Usadel equation Q 
with self consistency equation Q), we consider now this 
problem in the case of narrow domain walls in Sec. IIIII 
and large domain walls in Sec. IIVI 



III. NARROW DOMAIN WALL 

In this Section, we consider the case of thin domain 
walls characterized by the domain wall thickness w <C £ s . 
In Ref. 0, the zero temperature critical field was obtained 
in the context of magnetic superconductors. Here, we 
revise the result and obtain the phase diagram at finite 
temperature. 

We model the exchange field h acting on the electrons 
with a step function: h z (x) = /isgn(cc), h y = h z = 0. The 
structure of Usadel equation in spin space simplifies 
greatly and we have: 

- § dlf 11 + (M + is u h z (x)) f n = A, (10) 

while J 12 = / 21 = and / 22 = — /i 1 ,. Its solution for a 
given A is 



f n (x) = / dyg(x,y)A(y), 



(11) 



where Q is the Green's function associated with the ho- 
mogeneous differential equation 11U|) ; Q is defined by: 



G(x,y) = 
G(x,y) = 



Dk 



for x > y > 0, 



D(k + «*) 



K + K* 

e -Kx e n'y for x . > > (12) 



where k = + is u h)/D, while G{x,y) = G(y,x) 

and Q(x, y) = G(-x, —y)*- 

We look for a symmetric solution A(— x) = A(x). 
Writing the self consistency equation (lllll in Fourier 
space, we get the equation defining the critical tempera- 
ture T = T cw for DWS formation: 



(lnJL + 27rTRe^- 



1 



u> + ih + 



A„ = 



(13) 



u»0 J 



dk 



hVD[(u 



■h 2 ){Vu r 



h 2 



[h 2 + (w + ^) 2 ][^ 2 + (w + ^f) 2 



■A* 



Close to T c o, at h <^ T c q, the critical temperature for 
the transition into uniform superconducting state can be 
obtained analytically from Eq. ©: 



T cQ - T c (h) _ 7C(3) h 2 31C(5) h 4 



4^ 2 T 2 



167T4 T c 4 



(14) 



On the other hand, Eq. H3J) for the DWS can be simpli- 
fied: 



T cw (h) - T c {h) 7r Dp 2 



T, 



8 T, 



if) 



1 c0 



D 



<:■() 



dkA k , 
(15) 

where A = (8y/2— l)C(|)/(87r 3 ). Such equation is solved 
straightforwardly and we get the increase of critical tem- 
perature near T c0 due to DWS: 



T cw (h) - T c (h) 



1 c0 



(16) 



The corresponding shape of the order parameter near the 
transition is given by 



A(x) ~ cxp 



-B 



\x\ T cQ - T cw {h) 
6 



(17) 



where B = 167r v / 2A/(7C(3)). Thus, near T c0 , the lo- 
calized superconductivity is characterized by exponential 
decay without oscillation of the superconducting order 
parameter, with its maximum at the domain wall posi- 
tion. 

Away from T c q, Eq. I|13|) does not contain any small 
parameter. Its structure is that of a linear integral equa- 
tion whose kernel is a superposition of separable terms. 
Such form of the kernel is known to be convenient for 
numerical calculation. As a result, we obtained the sec- 
ond order critical line at any temperature (cf. Fig. 
The critical line is significantly increased compared to 
the critical line for the transition into uniform supercon- 
ducting state. In particular, at zero temperature, the 
critical field for localized superconductivity at T = 
is h cw — 1.33T c o ~ 0.76A and lies above the criti- 
cal field for the first order transition into uniform state, 
h c — 0.71Ao. It is of interest to note that, if the effec- 
tive exchange field in the bilayer is between h c and h cw , 
then the special situation occurs when only DWS can be 
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FIG. 1: Phase diagram. The critical line h^{T) at the second 
order transition into uniform superconducting state is plotted 
in black. At T < T* = 0.56T c o, transition into uniform state 
is of the first order, with the critical line h c (T) plotted in 
blue. The critical line h cw (T) corresponding to domain wall 
superconductivity is of the second order and is plotted in red. 
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FIG. 2: Self consistent gap just below the second order tran- 
sition line h C w{T) into localized domain wall superconduct- 
ing state at different temperatures T = 0, 0.2, 0.4, 0.6, 0.8 T c o 
(from narrowest to widest). 



realized in the system, but no superconductivity faraway 
from the domain walls. 

We also plot the self consistent order parameter at dif- 
ferent temperatures in Fig. [21 In addition to the decay, it 
shows small oscillations along the direction perpendicular 
to the domain wall at low temperatures. 

As it was pointed out in Sec. [HJ transition into uni- 
form superconducting state happens to be of the first 
order at T < T* in the bilayer. Therefore, one should 
also worry about the possibility for the change of the 
transition order along the critical line corresponding to 
DWS. We considered this possibility by solving the non- 
linear Usadel equations perturbatively up to the third 
order terms in the gap A. We found that such equation 
corresponds to the saddle-point of the free energy density 
(per unit thickness of the bilayer and per unit length of 
the magnetic domain wall) functional: 



T 



(18) 




UJ>0 • 



D 



11\2 



(/") 



where v is the density of states at the Fermi level in 
the normal state, and f^ 1 = (/I^)* where f* 1 is given 
by Eq. I|ll|) . At the second order transition, the term 
J- 2, which is quadratic in A, vanishes when Eq. I|13fl is 
satisfied. The phase transition is stable provided that the 
term J-^, which is quartic in A, remains positive along 
the transition line. We checked numerically that this was 
indeed the case. In particular, at T = and h close to 
h cw , we found: 



z/ln ■ 
h 

2.87v 



T A = 0.26z/ 



h 



dx\A(x) 



D 

2 fa/fti 



In- — A(x 

thr/tn 



0) 2 



D A{x = Qf 



2h r 



h 2 



(19) 



Thus wc found that the transition into DWS remains 
of the second order at all temperatures. 



IV. LARGE DOMAIN WALL 

In this section, we consider the case of a large domain 
wall with thickness w 3> £»• The domain wall can be 
described with an exchange field h = h(coscj), sine/), 0), 
the rotation angle cj>(x) varies monotonously between 
</)(— oo) = — 7r/2 and </>(oo) = ir/2. We find that the 
critical temperature at DWS nucleation is given by a 
Schrodingcr equation for a particle in the presence of a 
potential well whose profile is proportional to —(d x (f>) 2 . 
This result is not specific to the thin bilayer considered 
in this Section. In the Appendix, we extend it to bilay- 
ers with arbitrary thickness of the F layer and arbitrary 
transparency of the S/F interface. 

At the transition, the linearized Usadel equation Q 
must be solved, that is: 



-Y^/ 11 + H/ U -^(e-**/ ai - 
-?^/ 12 + M/ 12 + i/w w e- < */ 11 = 



-i<t> f!2\ _ ^ 



7 12 ) 



D 



9 2 / 21 + M/ 21 -i/ lS ^/ 11 = 0. 



(20) 



while f 22 = -f u . 

If the spatial dependence of <f> is neglected, the gap A 
is uniform along the layer and the solutions are readily 
found: 



u 2 + h 2 



-2i(f> f21 



pi 
JO i 



(21) 
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As a result, the critical temperature T c into uniform su- 
perconducting state naturally does not depend on (f> and 
is again given by Eq. ©. 

Now, let us assume that </> varies slowly. We solve 
Usadel equation (|20|l perturbatively by looking for a so- 
lution / f=a /o + fx- In first approximation, /o is still 
given by Eq. (|2T|) . where (f>, and, possibly, A, now slowly 
depend on x. The correction fx induced by their spatial 
dependence is determined by the set of equations: 

MA" - (e^fi 21 e-^fD = |^ 2 /o n 



Wl 2 +ihs„e-*f? = ^dlfl* 



(22) 



By appropriate linear combination of these equations, 
one finds that 



/i 2 



2(w 2 + /i 2 ) 5 



A' 2 A+ — — r^-rrA". 



2 (^ 2 + /i 2 ) 2 



(23) 

Here, primes stand for derivative along x. Inserting this 
solution in the self consistency equation Q we obtain 
the equation for the gap 

- — A"(x) + U(x)A(x) = EA(x), (24) 
2m 



where 



E 



In 



T 



1 \- 

— = dttt y 



(a; 2 - h 2 ) 
T P ' 2m ~ ^ (Z7TW : 



f/(x) = -DnT(rh') 2 ]T 



w>0 



(u 2 + /i 2 ) 2 ' 



(25) 



Equation (|24J) is a linearized Ginzburg-Landau equa- 
tion for a magnetic superconductor in the presence of 
a domain wall. It can be easily checked that effective 
mass m is always positive and U (x) is negative. There- 
fore, Eq. (|24|) looks like a one-dimensional Schrodingcr 
equation for a particle in potential well U(x). It is well 
known that a bound state with E < always forms in 
such potential. As a result, second order transition into 
localized superconducting state always appears more fa- 
vorable than into uniform state. 

Let us estimate now the magnitude of critical temper- 
ature increase. 

Close to T c o, the spatial variation for A is set by 
the temperature-dependent coherence length £(T) = 
£ s -\/T c o/ (T c o — T) which diverges at the transition. 
Therefore, £ s <C w <C £,(T) and the potential well 
can be approximated by a delta-potential: U (x) m 
-(7rL>/i 2 /48T c 3 w)(5(x), while l/2m - (-kD /9>T cQ ). There- 
fore we get the estimate: 



T — T 

-L CW J- C 

Tc. 



( h 



36 V27rT t 



cO 



e s e s iT cQ -T c f 



1 c0 



(26) 



Such increase is small in the ratio (£ s /ui) 2 , and another 
reduction factor comes from the smallness of critical ex- 
change field h near T c q, see Eq. {H}. 

At T — > 0, the second order transition into uniform 

(2) 

state occurs at exchange field h c = Ao/2. On the 
other hand, the effective mass in Eq. (|24|) diverges (as 
l/2m ~ (ttD /T)e~ h / T ). The fact that m remains pos- 
itive at finite temperatures is related to the absence of 
instability toward a modulated superconducting (Fuldc- 
Ferrell-Larkin-Ovchinnikov) state in magnetic supercon- 
ductors in the presence of strong disorder, tT c o <C 1^ 
Due to its large inertia, the particle now resides in the 
minimum of the potential well (zero point fluctuations 
can be neglected): C/ m ; n = — (irD /IQhw 2 ). The corre- 
sponding increase in critical exchange field at T = is 
(h cw - h ( c ] ) = itD/16w 2 . Such increase is of the order of 
magnitude ~ (£ s / w) 2 he <c h^. Actually, at low tem- 
perature the transition into the superconducting state in 
the uniform exchange field is a first order. We may expect 
that the domain wall superconductivity in such situation 
also appears by a first order transition. 

At intermediate temperature, we may obtain the crit- 
ical temperature T = T cw with a specific choice of the 
spatial dependence of rotation angle <\>. Assuming that 
4>{x) = 2 arctan[tanh(£/2tt))], we gettl*t 



In 



T 



8T r w 2 { 2irT r 



), 



(27) 



where 



F(u) =m,ij)x{Z) 




2$t(iip(Z)+wipx(Z)) 



(28) 

where i\) and t/>i are digamma functions, and Z = 1/2+iu. 

We should note however that transition into uniform 
superconducting state becomes of the first order at T < 
T* and results in significant increase of the critical line 
h c (T). Most probably, such change of the transition or- 
der should also be considered for localized superconduc- 
tivity. 



V. DISCUSSION 

The qualitative picture of the effect of domains walls 
on the superconducting properties now emerges from our 
calculations. 

First, the electrons of the Cooper pairs which travel 
across the S/F interface feel the exchange field hf in the 
F layer. This proximity effect results in an effective pair 
breaking that weakens superconductivity in S layer. As a 
result, the critical temperature T c of the bilayer gets sup- 
pressed in comparison with the critical temperature T c q 
of the bare S film: T c q — T c ~ i~ s _1 - There, the pair- 
breaking time r s can be estimated from Eq. (|14|) : at 
h <C T c q, t" 1 ~ h 2 /T C Q, where h is the effective exchange 
field that would act directly in the bare S layer to yield 
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the same pair breaking effect due to hf in the bilayer. In 
particular, when the S/F interface is transparent and F 
layer is thin, we found that h ~ hfdf/{d s + df), where 
df and d s are the thicknesses of F and S layers, respec- 
tively. On the other hand, when h S> T c q, there is no 
superconducting transition in the bilayer. 

Second, in the vicinity of magnetic domain walls in the 
F layer, such pair breaking mechanism becomes less effec- 
tive. Therefore, localized superconductivity may appear 
with critical temperature T cw > T c . When the domain 
wall width w is large in comparison with supcrcondutcing 
coherence length £ s , the exchange field rotates by the an- 
gle 8 ~ £ s /w on the scale of proximity effect. Therefore, 
the decrease of the average exchange field close to the wall 
is estimated as ft. — h av ~ 6 2 h. Correspondingly, the pair- 
breaking time is increased by At s /t s ~ 2 . In analogy 
with the theory of superconductivity at twinning-planc 
boundaries^, one can estimate the increase of T c : 

6 2 w 

AT = T cw -T c -7—r ■ (29) 

T s f(T) V ' 

Here, the temperature dependent correlation length 
£(T) ~ is \/T c /AT is the spatial extension of the su- 
perconducting gap and it diverges close to the supercon- 
ducting transition; the pair breaking is only reduced on 
the small portion of the gap corresponding to the width 
w of the domain wall. On the end, the formula yields 
the estimate AT/T C ~ (£ s / w) 2 / \t s T c ) 2 . Combining this 
result with the estimates for t s given in the preceding 
paragraph, we finally retrieve Eq. (|26[) qualitatively. This 
result holds when the width w is much larger than £ s . At 
smaller domain wall width, weakening of pair breaking ef- 
fect works on the characteristic scale £ s of the proximity 
effect and w should be replaced by £ s in the estimate. 
Therefore, the large enhancement of T cw , of the order of 
T c ), can be expected when w < £ s and h ~ T c0 . 

In the present work, we confirmed quantitatively this 
estimate in a number of situations. In particular, we 
obtained that, in the case of strong enhancement of T cw 
(at w < £ s ), for a thin F layer, transition into DWS 
state remains of the second order, with critical line above 
the transition into uniform superconducting state of the 
second order at T > T* = 0.56T c0 , and of the first order 
at T < T*. We predicted that only DWS could appear 
in bilayers with appropriate parameters. 

Let us now come back to the assumption of very thin 
F layer that was taken in the calculations. Usually, the 
exchange field in a fcrromagnet is much larger than the 
gap in a superconductor. Therefore, the coherence length 
£f is much smaller than £ s (1 — 5 nm for the former, com- 
pared to 10 — 50 nm for the latter). Therefore, the regime 
df <C £/ is quite hard to achieve with real samples of dif- 
fusive ferromagnets. Extending the calculation for thin 
F layer to arbitrary thickness of the layer and arbitrary 
transparency (characterized by 75) of the interface is 
quite straightforward when domain walls are large, as we 
show in the Appendix. As it is well known, the behavior 



of the critical temperature of the transition into uniform 
superconducting state is quite rich in such case and may 
even oscillate as a function of the parameters such as df 
or 7s. However, the physics of DWS appears to be quite 
similar to the one derived for thin bilayer. In particular, 
in the case of thick F layer (df £/ = y/Df/hf), the 
effective exchange field that would enter the above qual- 
itative estimates would be h ~ (£s/d s )y/hfT C Q. Clearly, 
at £ s ^> d s , as we assumed from the beginning, such field 
is much larger than T c0 and leads therefore to supercon- 
ductivity suppression. However, h is strongly reduced if 
the S/F interface is opaque, which may lead to the F/S 
coexistence and to DWS appearance, as studied in this 
work. Nevertheless, the enhancement of T c still is small 
by the factor £, s /w -C 1 in the situation described in the 
Appendix. 

These considerations suggest two possible directions 
to extend the range of existence of DWS. In the present 
work, DWS was analyzed for S layers with thickness d s 
much smaller than coherence length £ s . On the other 
hand, DWS should not appear when the superconduc- 
tor is hardly affected by proximity effect, at d s £ s . 
It would be of interest to consider the intermediate case 
when d s and £ s are of the same order. This was stud- 
ied for instance in the absence of magnetic domains in 
Ref. 0. Maybe a more important point would be to 
address the case of DWS in S/F bilayer with narrow do- 
main wall and large thickness of the F layer. However, 
both problems require considerably more numerical work, 
which goes beyond the scope of this paper. 

A lot of attention has been devoted to the study of 
long range triplet proximity effect which develops in S/F 
structures when the direction of exchange field in the fcr- 
romagnet varies spatially^SiiS. In our calculation, such 
long range triplet component is also present, as it is clear 
from the Appendix: in Eq. i|A6() . the term jochqoz in the 
direction transverse to the bilayer is generated only be- 
cause of the presence of domain wall, and it decays with 
typical length £r ~ y/Df/T £/. However, 70 does not 
enter f 11 and, therefore, is not important for the determi- 
nation of the critical temperature of transition into DWS. 
In all the calculations we presented, there is no long range 
triplet component in the direction transverse to the wall 
either: the typical length for the superconducting gap is 
determined by conventional short range proximity effect 
with decay length ~ £ s . We would like to emphasize also 
that the T c enhancement due to the appearance of DWS 
is maximized for narrow domain wall (see Sec. Illlfl . when 
the matrix elements of the anomalous Green's function 
that would give rise to long range triplet component are 
exactly zero. Therefore, the physics of DWS discussed 
here is not directly related to such phenomenon. 

This observation is in agreement with Ref. where 
the calculation of critical temperature of S/F bilayers in 
the presence of spiral magnetic order in the F layer was 
presented. There also, long range triplet component was 
found not to contribute to the result. On the other hand, 
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long range triplet component may be important for other 
properties such as density of state in the ferromagnetic 
layer with domain structure deposited on top of a bulk 
superconducting substrate^ 

We would like also to emphasize that our calculations 
differ from the study of S/F bilayers in the presence of 
spiral magnetic order in the F layer ^£ These works can be 
interpreted as considerations on magnetic domain struc- 
ture only in as so much that the width of the domains 
L and the width of the walls w are identical. A conse- 
quence is that the superconducting gap is spatially uni- 
form along the bilayer in Ref. fTst In contrast, our theory 
really shows that, in the more realistic case when w -C L, 
truly localized superconducting states can appear. In ad- 
dition, consideration on the effect of spiral magnetic or- 
der corresponding to 4>{x) = Qx in Eqs. I|24I25|1 . where Q 
is the wavevector of the spiral, can be immediately cal- 
culated from the Schrodinger-likc equation lj24|) . at least 
when Q^ s <C 1. Critical temperature enhancement due 
to spiral magnetic order and corresponding to uniform 
gap A(a;) follows straightforwardly from the observation 
that 4>'(x) = Q which enters Eq. (|25J) is constant. 

In conclusion, we analyzed in this work the enhance- 
ment of the superconducting critical temperature of su- 
perconducting/ferromagnetic bilayers due to the appear- 
ance of localized superconducting states in the vicinity 
of magnetic domain walls in the ferromagnet. We con- 
sidered the case when the main mechanism of the super- 
conductivity destruction via the proximity effect is the 
exchange field. We demonstrated that the influence of 
the domain walls on the superconducting properties of S 
layer may be quite strong if the domain wall thickness is 
of the order of superconducting coherence length. 

We interpreted qualitatively and quantitatively the 
amplitude of this effect, and we pointed out the special 
case when parameters of the bilayer arc such that only 
localized superconductivity may form in these systems. 

For magnetic film with perpendicular anisotropy, the 
orbital effect provides an additional mechanism for the 
domain wall superconductivity* and it may be easily 
taken into account. On the other hand for the film with 
the easy plane magnetic anisotropy the domain wall will 
be a source of the magnetic field in the adjacent S layer 
and locally weakens the superconductivity^ The role of 
the orbital mechanism in the domain wall superconduc- 
tivity may be important only if the magnetic induction 
is comparable with the upper critical field of the super- 
conducting film. 

The domain wall superconductivity in S/F bilayers 
opens an interesting way to manipulate the supercon- 
ducting properties through the domain structure. In par- 
ticular the motion of the domain wall in F layer may be 
accompanied by the displacement of the narrow super- 
conducting region in S layer. 

We are grateful to J. Aarts for attracting our attention 
to the problem of the domain wall superconductivity in 
S/F bilayers and useful discussions. 



APPENDIX A: THICK F LAYER 

As mentioned previously, hf is usually large in fer- 
romagncts. Therefore the condition of thin F layer, 
df <C £/, is hardly reached. We would like to extend 
the results of the previous section to the more realistic 
situation of a finite size F layer. The difficulty is that the 
set of differential equations Q-lQJ to be solved are now 
two-dimensional. We managed to solve it for the case of 
large domain wall only. In the F layer, we parametrize 
the exchange field rotation with a slowly varying angle 
<j>(x) such as hf = hf (cos </>, sin cj>, 0) . Linearized Usadel 
equations @ in F layer are: 

A/ 11 + M/ 11 - (e-^/ 21 - e"^/ 12 ) = 0, 

-^A/^ + M/^ + ifey^e-*/ 11 ^, 



D 



f a -p21 



Af ZL + \iu\f ZL -ihf S ^f lL =0. 



(Al) 



where A = <9 2 + <9 2 , while / 22 = — / n . When the spatial 
dependence of 4> is neglected, the general form of the solu- 
tions which satisfy boundary condition at the F/vacuum 
interface are: 



12 



/o 1] 
/o 
/o 21 



F + chq + z + F_chg_z, 



(A2) 



s u e l<t> {Fochgo-S + F + chq + z — F_chq_z} , 
—Sue {— -Fochf/oz + F + chq+z — F_ch<7_z} . 



where z = z + df, qo = ^/2\uj\/Df, and q± = 

y/2(±ihf + \u\)/D f . 

We determine now the amplitudes of the cigenmodes 
Fq and F±. For this, we first determine f 11 and <9 z / n at 
z = in F layer. Making use of the boundary conditions 
(0), wc can now insert them in the Usadel equation © 
in the S layer. We proceed similarly for f 12 and / 21 . On 
the end, we get from JSJ) three equations which determine 
the amplitudes we are looking for. We find Fq = 0, while 



F± = 



A 

2TU 



fl± = \u)\C± + aq±shq±df, (A3) 



where C± = chq±df + "fB£,s<l±shq±df and a = 
D s af /2d s a s . 

Inserting Eqs. (|A2(I - I|A3(I into J2J to determine / s , and 
then inserting / s n in the self consistency equation Q), we 
get the equation defining the critical temperature T = 
T c (h) for uniform superconducting state: 







ln^+27rTRe^ 



r+ = 



aq + 



coth q+df +7 B £ S <7 + 



(A4) 



This result is the same as Eq. (47) of Ref. fl Whether 
transition is of the second order (as described by 
Eq. (|A4p ) or of the first order was considered in Ref. |2j 
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Let us now consider the effect of a domain wall. As in 
section II VI we will look for a solution / w f + fx . In 
leading order, the spatial dependence of <j> and the gap 
A is ignored, and f is still given by Eqs. (|A2IA3|I . The 
correction f\ accounts for the slow variations of <f> and A 
along x-axis; it is determined by the set of equations: 



Finally, we can insert Eqs. i|A7(l into HA6jl . and then 
into in order to obtain f 11 in the S alyer. Then, we 
insert it the self consistency equation On the end, 
we find that the superconducting gap at the transition 
into DWS state is still determined by the Schrodingcr 
equation 1(23}, with effective coefficients: 



^dlfl 1 + 1« 


l/r 11 -* 


hfS^ , _ 
2 (e 


-i(j> f 21 
Jl 


(A5) 






— e~ 




- ^d 2 f n 


-^.A 12 + 


MA 12 h 


- ihfS^e 


— i<p rll 
Jl 


- ^d 2 f 12 


-^A 21 - 


f- MA 21 


- ihfS^ 


e i<f, f n 


- 2 °*J° ■ 



Ignoring ^-dependence of the right-hand-side of the 
above differential equations, we can find the exact in- 
dependent function f\ which solves them. We obtain: 



E 

«=± 



A iD f( j)' 2 A 
chq n z — — — h a — - — — — 
1 2Ci a WhfCl a 



la 



-zshq a i 



A" 



\<f>>*A 



a=± 



iD f 



16hfQ a 



'■A D f (2<f/A' 



aA 

2n a 

'A) 



-azshq a z- 



A'' 



4:h f fl a 
A" 



+ aia 
c\iq a z 
(A6) 



where the integration constants 70 and r y± still remain 
to be determined. For this purpose, we insert Eq. (|A6|) 
into (JH} in order to get f s and d z f s . Inserting them in 
Eq. JS}, we thus obtain a set of three equations which 
allow to determine them. In particular, we find: 



7± 



± 



iDf(j)' 2 A 
16h f n ± 

4 U 2 ± T 



" ± (A" 



/ 2 A) 



i<j>' 2 AC ± 



7^ (A7) 



where E± = |oj|*S'-i- +a(shg±iij + dfq±chq±df) and S± 
dfshq±df + 7sCs(shg±d/ + d/g ± chg ± c?/). 



£ = 



(A8) 



1 

2m 



^srV^ , Oi(sh2q + d f + 2q + df) 



UJ>0 



U(x) = -ttT(0') 2 E ^ 



w>0 



„ 1 

0J+ 



2q+n 2 + 
2 



where w 



+ 



3? a(sh2g + d / + 2 g+ d/) _ _D£^J_ 



In the limit of thin F layer {q+df — > 0) and large 
transparency of S/F interface (75 — ► 0), we note that 
r + ~ aq 2 + df ~ (rjf/rj s )hf, and it is easily checked that 
the above formulas for l/2m and U(x) reduce to Eq. I|25(l 
from Scc.^l when rjf <C r] s . 

For large F films and transparent S/F interface, we 
find that the critical temperature is given by Eq. l|Al} . 
where T + = (1 + i)h and h = a.\fhf / Dj: T + can be in- 
terpreted as a combination of both exchange field and 
spin-flip terms with equal weight. In such case, it is 
known that transition into uniform state is of the sec- 
ond orderiSi The coefficients of Eq. l|A8|l also simplify to 
the form 



lr 2m 



1 



w>0 



[u+ (1 + i)h]< 



U(x) = -nTD^'fY 



UJ>0 



1 



uj+ (l+i)h 



(A9) 



Let us note that effective exchange field and spin-flip pa- 
rameter scale is h ~ {^ s /d s )^JhfT^a if S and F layers 
have comparable diffusion constants and conductivities. 
When d s <C £ s , as we assumed from beginning, this leads 
to h~^> T c q , and therefore to complete superconductivity 
suppression. However, we expect that our results hold 
qualitatively in more general case £ s < d s when super- 
conductivity is not completely suppressed. 
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